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fields with a non-trivial moduli space metric. We discover that the scaling property of the 
cosmology is synonymous with the scalar fields tracing out a particular class of geodesies in 
moduli space - those which are constructed as integral curves of the gradient of the log of the 
potential. Given a generic scalar potential we explicitly construct a moduli metric that allows 
scaling solutions, and we show the converse - how one can construct a potential that allows 
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1 Introduction 

Scaling solutions to the Friedman equation are of prime interest in modern cosmology, both for 
study of early time inflationary scenarios and, more recently, the late time accelerating expansion 
of the universe which is suspected to be driven by some dark energy component in the cosmological 
fluid 1^ 121 . The utility of these scaling regimes is that the various components can evolve such 
that a constant ratio is maintained between their energy densities, and as attractor solutions they 
approach these ratios without fine tuning of the initial conditions and so can be used to explain the 
cosmic coincidence of energy densities 011]. Scaling solutions also find a use in the earlier epoch of 
infiation, where exponential potentials |S1 El Cj allow for exact solutions and sums of exponentials 
can aid infiationary behaviour 8 . The dynamics of scalar fields with such potentials is now well 
understood, and is neatly expressed using ideas from dynamical systems for single scalars (HI and 
multiple scalars [inilTTlIT^. 

That scalar field models can be constructed with potentials that can mimic the conditions 
we observe in the universe is but one part of the story. The theories proposed require a fine- 
tuning to match observation, and we are still as of yet to detect the direct signature of scalar 
fields in nature; is there any fundamental physics that might motivate their existence? Fortunately 
scalar fields abound in the low energy limit of many modern particle physics models where they 
appear dynamically in the form of moduli which parametrize the size and shape of the geometry 
of internal dimensions. Interestingly in many of these models, particularly unified gravity models 
with dimensional reduction, gaugino condensation or instanton corrections, the low energy effective 
potentials are often exponential in form T31|Ill|Tni|Tni|I71|THllIlllinillTl|121- It would be intriguing 
indeed if the scalar potentials used in infiationary theory and quintessence were in fact a consequence 
of dynamics of some extra dimensional string or M-theory. 

However, exponential potentials are but one class of possible potential. How might we inves- 
tigate the scaling behaviour of more general potentials? As cosmologists we are well aware of the 
geometric nature of pure gravity in Einstein's equations, however when we couple in additional 
scalar fields they also come with their own geometry. As the system evolves one can think of the 
kinetic term providing a measure of distance in field space. This is by virtue of the existence of a 
moduli space metric, which for minimal kinetic terms is flat, d^(f)"'d^(j)'^ = dabd^(/)"'d'^(j) , but there 
is no reason a priori not to also consider other forms for this metric too, perhaps providing a useful 
geometric insight into the dynamics of the system. 

There is already a tradition of the study of moduli space metrics within particle physics models, 
particularly in supersymmetric held theories where the field space geometry is usually of some 
special type, for example Kahler geometry. More recently Townsend and Wohlfarth proposed 
thinking of cosmological evolution in an augmented field space |231 1 where the scale factor is added 
to the geometry of the field space. It is then discovered that cosmological solutions can be thought 
of as geodesies on this extended manifold. Here however we do not consider this augmented space, 
rather we look at the field moduli space alone and discover that in a scaling regime the fields 
trace out a particular type of geodesic, namely a geodesic which coincides with an integral curve 
of daln(y). Given that the moduli metric. Gab, and the potential, V, are in principle independent 
there may be no such geodesic, or indeed there may be many congruences of such geodesies. Here 
we explicitly construct a Gab (although there may also be others) that allows scaling solutions with 
a generic potential and we also show how one can find a potential that allows scaling if Gab is 
known. 



2 Scaling Cosmology and Geodesies 

The starting point for this work is a simple observation following from previous studies of multi-field 
systems with Gab = ^ab- It has been shown ^l^^E] that the scaling regime is characterized by 

r = — , (2.1) 

^0" = A''ln{t) + B''. (2.2) 

So given two different fields we have 

which is just the equation of a straight line. As the moduli metric is flat for these cases then this 
straight line is, of course, a geodesic. It is our aim to show that even for Gab ¥" ^ab this relation 
between scaling and geodesies persists. 

We start by describing the system we aim to study, namely a set of scalar fields with moduli 
metric Gab and potential V coupled to gravity. The Lagrangian is given by 



C = V^g 



^n-\Gabd^4>''d^^'' -y{^) 



(2.4) 



with the indices a,b,c,... G {l,2,...,n} running over the number of scalar fields. We adopt a 
conventional Friedman-Robertson- Walker cosmology^ for the space-time metric and restrict the 
scalar to vary only with time. This then allows us to write down the effective Lagrangian for the 
reduced system, where the dynamical variables depend only on time, 

C = a' Qg,50>' - 1/(0)) - ^ah\ (2.5) 

Here a = a(t) is the scale factor of the universe, and ' = dt = ^ denotes the derivative with respect 
to time. In this language we find that the equations of motion, along with the vanishing of the 
Hamiltonian due to invariance under time reparametrization, are 

H^ = '^(^^Gabr^' + Vy (2.6) 

H + 3H'^ = kV, (2.7) 

4>'' + T%J^<P^ + 3H<j>'' = -G''''dbV, (2.8) 

with the Hubble parameter H = - and the affine connection r^,^ of Gab defined in the conventional 
way, 

r^c = ■^G"^ {Gbd,c + Gcd,b — Gbc,d) ■ (2.9) 

Given that we shall be discussing geodesies we now give the standard expressions for a geodesic 
curve, C, with tangent vector T that is parametrized by an affine parameter A. 

r« = -^</>^ (2.10) 

dX 
T^VaT^ = ^ (p"'' + T%^(p'^(P"' = 0, (2.11) 

T^Ta = 1, (2.12) 



^In this paper we take a space-time metric of signature ( — \- ++). 



where ' = ^. 

We shall take scaling to mean that the terms in ()2.6() and (|2.7|) evolve in constant proportion 
to one another, H^ (x H (xV (x Gab4''^4' ■ We start by considering H^ and the kinetic term of (|2.H|) 
and enforce proportionality by introducing a constant a, 

a^H^ = Gabr^'. (2.13) 

We now want to see if this can be consistent with (p^{t) tracing out a geodesic, meaning that the 
affine parameter can now be thought of as a function of cosmic time A = A(t). The normalisation 
of the tangent vector (|2.12|) then gives us 

Gabr^'=(^) , (2.14) 



\dt 
which combines with ()2.13|) to produce 

a = oo exp(A/a). (2.15) 



Converting time derivatives into A derivatives and using (|2.7)) along with (|2.1H) one finds that (|27 
becomes 



^ C^dbHV). (2.16) 



ok"^ 



This is very suggestive, and to make manifest its relation to geodesies more explicit we shall write 
it as the components of the tangent co-vector, 

Ta = daW, W = ^ln(y). (2.17) 

So, if our scaling/geodesic ansatz is consistent with the other terms in the equations we see that 
the scaling geodesic takes a very particular form, it is the integral curve of VW. An immediate 
consequence of this is that one can never get a scaling solution on closed geodesies. 

We also need to check that the potential evolves in constant proportion to H^, which we can 
do by using the fact that these tangent vectors have been normalized to unity, 

(t>"'{-^)daHV) = l, (2.18) 



an"^' 



which is integrated to give 



V = Voexp{-aK^X) = Vo[ — ] . (2.19) 



If it is the case that V ex H^ then we can integrate to find 

a"'^y^oc{t-to), (2.20) 

which, when substituted back into the potential ()2.19l) . gives the conclusion that 



1 



yoc-2. (2.21) 



This is precisely the form that V must take as is easily seen by inspecting ()2.6()(|2.7() . We also 
note from ()2.2U|) that we get an accelerating scale factor if ok < \pl. We can understand such a 
requirement by reconsidering (|2.12|1 and (|2.17|) and writing them in the form 

^G^^'daVd^YlV'' = \a\\ (2.22) 

and now we recognize the left hand side to be the multi- field generalisation of the slow-roll parameter 
e 121] • We should point out that the potential is only required to take the form H2.19() on the geodesic 
itself, away from the geodesic the potential can take any form, or even be undefined. 
A brief comment is in order about the nature of a. We see from ()2.6p and ()2.13() that 

and so it seems in scaling solutions a determines the value of both the kinetic and potential energies. 
Moreover, if the potential is positive then we have an upper bound on a given by aK < v6- 

3 Construction of a field space metric 

We now come to the issue of whether we can construct a metric Gab that will allow scaling for a 
given potential V. Using 1)2. 17() along with 1)2. 11() we find that this amounts to solving 

C'^w (dbdcW -ri.ddw'^ = o, (3.25) 

for the metric. At first sight this looks a difficult task, and a general solution probably is out of 
reach. However, it is possible to solve this equation if we restrict ourselves to a conformally flat 
metric 



ds^ = exp(2n)d^% (3.26) 



with connection components 



r'fec = S^^,c + 5^^,b-5bc6"' ^,d- (3.27) 



This is a natural choice in that (|3.25|) is now an equation for a single scalar function 0,{(p) in terms 
of the single scalar function V{(f)). Given that we have already normalized the tangent vector in 
(|2.12|) this immediately gives us the conformal factor 

exp(217) = Y. daWdaW = J^ J: ^^ (3.28) 

a a 

which, as is easily checked, does indeed satisfy (|3.2f)|) . One consequence of this solution is that if V 
goes through zero with non-vanishing gradient then the conformal factor diverges there, pushing 
such points to an inflnite proper distance in moduli space. Hence in a scaling solution V will 
not change sign. Another important aspect of this solution is that it allows for a scaling solution 
through each point in field space where ln(y) is well defined; it does not restrict to a single path 
as is the case for the usual multi-field scaling solution. 



The conformally-flat solution we have presented is not necessarily the only metric on field space 
to give scaling solutions, however we leave investigation of the construction of more general metrics 
for future work. 

The other way to proceed is to find a potential that yields a scaling solution for a given metric. 
We do not have an explicit solution to this problem but one certainly exists. Once a moduli metric 
has been chosen one finds a suitable, i.e. not closed, geodesic C. Using the affine parameter of 
this geodesic one can reconstruct the value of the potential along C using 1)2. 19() . Simply having V 
along the geodesic is not enough to keep (/>" on the scaling solution, we also need daV which we 
can find from (|2.17)) as we know the geodesic's tangent vector. This is rather like a race track with 
banked corners to keep the cars from coming off - the gradient of V keeps the scalars on the scaling 
solution. 

4 Constructing scaling solutions 

Now we have set up the framework for finding scaling solutions in general potentials we shall give 
some examples to show how the process works, by numerically evolving (|2.6j)(|2.7))(|2.8j) . We shall 
start with a rather simple example, reminiscent of the multi-field scaling solutions which required 
exponential potentials. 



F = exp(0i) + exp(2(/)2), 



(4.29) 



using the values k^ = 1 and a = 1. In Fig. ^we have plotted the evolution of (p^ and cf^, also 
showing the vector field da In V (the arrow lengths have been normalized for clarity) . The solid 
lines denote scaling solutions and we note that the fields do indeed follow the gradient of ln(y) as 
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Figure 1: This plot shows the vector field Vln^ of the potential (|4.29|) . along with some example 
field evolutions. The solid lines are the scaling solutions, with the dashed lines being perturbations 
about the scaling solution. 
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Figure 2: This plot shows the ratio of the kinetic, potential energy to H^ for a scaling solution 
and perturbations around it. 




Figure 3: This is the gradient of the potential (|4.H()|) . along with the scaling solutions. 



expected. The dashed Unes represent perturbations of these initial conditions, and we see that they 
do not simply go back to the scaling solutions of the solid lines. This however does not mean that 
scaling solutions are not attractors as Fig. [21 makes clear. Fig. |21is a plot showing the values of the 
kinetic and potential energies for a scaling solution (solid line) and a perturbation of that solution 
(dashed line). What we see is that the perturbed case does approach a scaling solution, just not 
the one you started with. This is because when using the conformally flat moduli metric to solve 
H3.25() there is a scaling solution through each point in moduli space. Note also that T/H'^, V/H"^ 
approach the values as predicted by 1)2. 23(1 and ()2.24|) . 

As a further test let us consider a more involved, bumpy, potential, 

V = (2 + sin(0i))(2 + sin((^2))^ (430) 

Scaling solutions for this can be seen in Fig. |21 also following integral curves of the gradient of 
ln(y), as our analysis shows they should. 

5 Adding a barotropic fluid 

The effects of adding a baratropic fluid to the system can be studied by extending our Friedman 
equation (|2.(ij) and adding the fluid equations, 

^2 ^1 



H' = —{-Gab(t>''(t>' + V + p\, (5.31) 

P = (7-l)p, (5.32) 

p = -3'yHp, (5.33) 

where P is the fluid pressure and p is its density. Such a fluid is typically invoked to represent 
species of particles in the early Universe other than the scalar field and so constitutes an important 
addition. 

Although we have not thoroughly explored the consequences of adding such a fluid to the system 
we can find a bound on the equation of state parameter, 7, which reduces to the standard one of 
the usual scaling analysis. The analysis follows the same lines as the one presented above, using 
the scaling ansatz H2.13|l as well as 

V + {l-j/2)p = pV, (5.34) 

and taking /3 as a new parameter which is constant during scaling. With this addition we find that 
1)2. 16() is altered to become 

^ C^dblniV). (5.35) 



a(3K 



A consequence of 1)5. 34() is that if the fluid has 7 < 2, which is physically reasonable, then for 
positive V and p we require (3 > 1. This is the origin of the bound we shall derive. 
In much the same way that (|2.19|1 is derived, one finds that 



p oc a-" ^^ , (5.36) 

however we also know by integrating ()5.33|) that p oc a^^'*' which gives us 

P = 3-f/{a^K^), (5.37) 



7 



and so our bound on (3 now becomes a bound on the equation of state parameter 

37 > a^K^. (5.38) 

This bound is similar to the one usually found in scaling analyses with a fluid, indeed taking the case 
of a single scalar with potential V = Vq ex.p{bK(/}) one has that o^k^ = 97^/6^, so 1)5. 38(1 becomes 
6^ > 37 matching the result found in jH]. 

6 Comments on fixed point analyses 

No paper on scaling dynamics appears to be complete without consideration of the critical point 
analysis, obtained by recasting the equations of motion into an autonomous systems framework 
[ni^J. In this approach the asymptotic behaviour and stability of the system can be discovered 
by analysing the nature of fixed point solutions. 

The now standard approach for a system with potential 

V = Y,Aiexp{-ai.^) = Y,Vi (6.39) 

is to define two new sets of variables 

X" ~ ^""/H, (6.40) 

Y^ ~ V^/H, (6.41) 

with which to describe the degrees of freedom. These are useful because in the case of exponential 
potentials they are constant in the scaling regime, giving scaling solutions as fixed points of the 
autonomous system. For the more general case we have been discussing the closest analogue would 
be the term 



X ~ ^^Gabr(t>'/H, (6.42) 

which we know is a constant during scaling. However, this is just one term and we need to describe 
n momenta (p"". In the case where one uses the conformally flat solution of ()3.25|) one could introduce 

X" ~ 6^4)" /H, (6.43) 

as a natural set of variables. However, the problem reappears when we try to find an analogue for 
the Yi. In general the potential will not split up nicely into a sum of terms, making the analysis 
problematic. One possible approach to constructing more variables is to use the derivative of the 
potential, this is a problem which we hope to return to. 

7 Concluding remarks 

In this paper we studied the scaling solutions of a Friedman- Robertson- Walker cosmology containing 
scalar fields evolving in an arbitrary potential. The standard picture with canonical kinetic terms 
would lead us to conclude that the potential must take a special form, that of a sum of exponentials, 
if scaling is to take place. What we have been able to show is that by including the freedom of a 
metric on moduli space we can recover scaling solutions in generic potentials, with the connection 
between moduli space geometry and scalar potential being that the scaling geodesies lie along 
integral curves of VlnT/. In fact we have been able to provide a special solution to the scaling 



criteria by way of a confornially flat moduli metric, and we provided explicit examples to show the 
system in action. The statement also works in reverse; given a metric on moduli space it is possible 
to construct a potential which will support scaling solutions. We also briefly considered the effects 
of adding a baratropic fluid into the system and were able to find a bound for the equation of state 
parameter, with this bound reducing to the standard one when the moduli metric is flat. 

There are a number of things which have not been fully addressed is this article. Firstly is the 
issue of stability. While we saw numerical evidence in our examples that the system approached a 
scaling solution a more thorough analysis is required, and we expect it will combine some aspects 
of the moduli geometry along with the potential. An approach that could help with the stability 
analysis is that of recasting the system in a form inspired by dynamical systems, writing the 
equations as an autonomous system. While we have not presented such a framework we believe 
this approach is well worth studying. 
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